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Notions of a Gaussian state and a Gaussian linear map are generalized to the case of anticom- 
muting (Grassmann) variables. Conditions under which a Gaussian map is trace preserving and 
(or) completely positive are formulated. For any Gaussian map an explicit formula relating cor- 
relation matrices of input and output states is presented. This formalism allows to develop the 
Lagrangian representation for fermionic linear optics (FLO). It covers both unitary operations and 
the single-mode projectors associated with FLO measurements. Using the Lagrangian representa- 
tion we reduce a classical simulation of FLO to a computation of Gaussian integrals over Grassmann 
variables. Explicit formulas describing evolution of a quantum state under FLO operations are put 
forward. 



PACS numbers: 03.67.-a 

I. INTRODUCTION 

Photons and electrons are fundamental particles which 
are widely used in the modern communication and com- 
puting technology Since these particles can be used 
to represent quantum information as well, they could 
be involved in quantum computing technologies of the 
future. It was discovered recently by Knill, Laflamme 
and Milburn p| that passive linear optical elements to- 
gether with squeezers, single photon sources and photo- 
detectors provide a universal set of operations for quan- 
tum computation. After that discovery Terhal, DiVin- 
cenzo 0, and Knill Q examined a computational ca- 
pability of "fermionic linear optics" (FLO). This term 
basically refers to a system consisting of non-interacting 
electrons in a controllable external potential and a de- 
tector that measures projectively occupation numbers of 
single-electron modes. It turned out that FLO is not 
a promising platform for building a universal quantum 
computer, since it can be efficiently simulated by clas- 
sical means. On the other hand, classical simulatibility 
makes FLO interesting from the perspective of quantum 
communication theory, since it might help to understand 
better general properties of quantum channels. 

As we shall argue, FLO shares many features with 
a restricted version of photon linear optics where only 
homodyne measurements are allowed. Using linear op- 
tical elements, squeezers, and homodyne measurements 
one can generate only Gaussian states, and realize only 
Gaussian maps, see 0, which also can be simulated 
classically. One reason for this similarity is that non- 
interacting fermions (bosons) are described in the La- 
grangian representation by Gaussian integrals over an- 
ticommuting (commuting) variables. In the present pa- 
per the notions of a Gaussian state and a Gaussian lin- 
ear map are generalized to the case of anticommuting 
(Grassmann) variables. Conditions under which a Gaus- 
sian map is trace preserving and (or) completely positive 
are explicitly formulated. 

As was pointed out in Q , the weakness of fermions is 



that all FLO operations, both the unitaries and the pro- 
jectors, belong to a closure of a Lie group of relatively 
small dimension (growing only quadratically with a num- 
ber of modes). We shall call elements of this Lie group 
as Gaussian operators. Roughly speaking, any Gaussian 
operator can be represented as an exponent of another 
operator which is quadratic in the creation/annihilation 
operators. Unfortunately, this definition can not be ap- 
plied directly in many important cases, when an oper- 
ator does not have a full rank. We show how to avoid 
problems related with taking a closure and considering 
limiting points by defining Gaussian operators in terms 
of anticommuting variables and taking advantage of their 
symmetry properties. It turns out that the adjoint action 
of any Gaussian operator is a Gaussian linear map. 

A set of states that can be achieved by FLO opera- 
tions starting from the Fock vacuum is a set of fermionic 
Gaussian states. As it is the case with the bosonic coun- 
terpart, a fermionic Gaussian state can be completely 
characterized by its correlation matrix and thus admits 
an efficient classical description. We shall study an ac- 
tion of Gaussian maps on Gaussian states in details and 
describe explicitly how correlation matrices of input and 
output states are related to each other. 

A classical simulation algorithm for FLO must be able 

• to reproduce statistics of all measurement out- 
comes, 

• to describe explicitly a quantum state of the system 
at each time step. 

The first issue has been addressed in the paper , where 
an explicit formula for the measurement outcomes proba- 
bility distribution has been put forward. In a more recent 
paper 6] DiVincenzo and Terhal described a transforma- 
tion of FLO states under single-mode measurements in 
a special case, when the total number of particles in the 
system is conserved. The representation of FLO in terms 
of Gaussian maps allows to describe in an explicit form 
a transformation of correlation matrices under the whole 
set of FLO operations. 
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The rest of the paper is organized as follows. Sec- 
tion [n] contains a strict definition of FLO. Section ITTT1 is 
a brief outline of the Grassmann variables formalism. It 
is based primarily on the textbook and the paper @- 
In Section IPTl we use Grassmann variables to introduce 
a convenient parametrization of linear operators and for- 
mally define fermionic Gaussian states. A set of Gaus- 
sian operators incorporating Gaussian states, canonical 
Bogolyubov transformations, and FLO projectors is de- 
fined in Section Theorem ^ therein provides a non- 
singular characterization of Gaussian operators. It proves 
that FLO transformations preserve the set of Gaussian 
states. A Lagrangian representation for linear maps (i.e. 
transformations sending operators to operators) is devel- 
oped in Section lv"fl where we define a notion of fermionic 
Gaussian map. Conditions under which a map is trace 
preserving and bistochastic are described. Completely 
positive Gaussian maps are characterized in Section fV 111 
Finally, in Section IVTTT1 we describe Gaussian maps im- 
plementing FLO operations and find their action on the 
level of correlation matrices. 

II. FERMIONIC LINEAR OPTICS 

We start from giving a more strict definition of FLO. 
Instead of talking about electrons, it is more convenient 
to define abstract fermionic modes which can be de- 
scribed by creation and annihilation operators a] and a,j , 
j = 1, . . . ,n, with n being the total number of modes. 
They satisfy the Fermi-Dirac commutation relations, i.e. 
{a,j,ak} = and {a.j,a\} — Sj t k- An occupation number 
Nj of the j-th mode becomes an observable Nj — a^dj. 
The vacuum state jO^) satisfies a,j\0p ) = for all j. The 
Fock basis is defined by 

\N u ...,N n ) = (a\) Nl ■■■(al) N " \0 F ). 

Here the occupation numbers take values Nj £ {0,1} 
according to the Pauli principle. An arbitrary state of 
the system can be written as a superposition of the Fock 
basis states. 

Quantum computation with FLO starts from a Fock 
basis state (w.l.o.g. from the vacuum state). It is acted 
on by a sequence of "linear optical" elements. Unitary 
elements describe an evolution under quadratic (in cre- 
ation/annihilation operators) Hamiltonians, that is 

n 

H = Y J ^a)a +H t + H s (1) 

where Cj represent individual energies of the modes, 
H t = 2J tjkOj&k +tjka' k aj 

l<j<k<n 

represents "tunneling" processes, and 

H s = ^2 Sjk&'j&l + Sjkakaj. 

l<j<k<n 



If the fermions under consideration are electrons, H s may 
describe an interaction between the system and a bulk p- 
wave superconductor In this case creation of two 
extra electrons in the system is balanced by borrowing 
one Cooper pair from the superconductor. 

Non-unitary "linear optical" elements are single-mode 
measurements. We can projectively measure an occupa- 
tion number Nj for any mode j. Depending upon the 
outcome, the input state is acted on by a projector OjCij 

or by a projector fijoj. In a FLO quantum computation 
we can alternate between measurements and unitary el- 
ements, such that a choice of the next operation may 
depend upon the prior measurement results. 

Instead of 2n creation and annihilation operators it 
will be convenient to introduce 2n Hermitian operators 

C2j-i = fit + &j and c 2j = (— i)(a] - a ). (2) 

They are analogous to coordinate and momentum oper- 
ators for bosonic modes. It follows from the Fermi-Dirac 
commutation relations that 

c a Cb + c b c a = 2S ab I (3) 

for all 1 < a, b < 2n. An algebra generated by the opera- 
tors {ca} is called the Clifford algebra and will be denoted 
C2n- An arbitrary operator X £ C211 can be represented 
as a polynomial in {c a }, namely 

2n 

X = aI + ^2 X] a ai,...,a p C ai ---C ap , 

p—1 l<ai<---<a p <2n 

where a = 2~ n Tr (X). We shall call an operator X £ €2,1 
even (odd) if it involves only even (odd) powers of the 
generators c. Thus the operator algebra describing n 
fermionic modes coincides with C2„. In particular the 
Hamiltonian in Eq. looks as 

2n 

H = ~ X H ab c a c b , (4) 
a,b=l 

where {H ab } may be an arbitrary antisymmetric real 2nx 
2n matrix. Unitary elements of FLO constitute a group 
of canonical transformations G c C U(2 n ), that is V £ G c 
iff V = exp (iH) with H from Eq. Adjoint action of a 
canonical transformation V £ G c on the operators {c a } 
is just a rotation: 

2n 

Vc a V^ = J2 R abC b , RR T = I, det(Ji) = l. (5) 

b=l 

By choosing appropriate H one can implement an arbi- 
trary rotation R £ SO(2n). 

III. ANTICOMMUTING VARIABLES 

It has been known since the late seventies that a pow- 
erful technique to study systems with infinite number 
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of fermionic modes (fermionic quantum fields) is the 
functional-integral, or Lagrangian representation (LR). 
In this representation many interesting quantities such 
as the Green's function or the partition function are ex- 
pressed in terms of Feynman integrals over Grassmann 
valued field. A concise exposition of this technique can 
be found in the original paper by Soper 8] . Of course LR 
can be applied to a system with finite number of modes 
as well. 

In the present paper we slightly customize LR to make 
it more suitable for description of FLO. First of all, since 
the unitary elements of FLO do not preserve the total 
number of particles, we do not have a preferred vacuum 
state. Accordingly, we do not have a preferred normal 
ordering of fermionic operators. Secondly, we would like 
to treat both pure and mixed states on equal footing. It 
makes more convenient to develop LR for mixed states 
and for linear maps sending mixed states into mixed 
states, rather than for pure states and for linear oper- 
ators sending pure states into pure states. For the sake 
of completeness, we start from a brief outline of the an- 
ticommuting variables formalism, see a textbook |7j for 
more details. 

Consider n-dimensional complex linear space C n and 
denote 9±, ...,6 n its basis vectors. Grassmann algebra 
with complex coefficients Q n is generated by 9\ , . . . , 9 n 
considered as formal variables subject to multiplication 
rules 



9^ = and 9 a 9 b + 9 b 9 a = Q. 



(6) 



An arbitrary element of / G Q n can be written as a poly- 
nomial of 6*'s, namely 



/(*) = «+£ E 



*«1 , . ...a p u a\ 



p—1 l<ai<---<a p <n 



where the coefficients a* are complex numbers. In other 
words Q n is an antisymmetric tensor algebra over n- 
dimensional complex space. A polynomial f(ff) will be 
called even if it involves only even powers of 9. Even el- 
ements constitute the center of the Grassmann algebra. 

One can formally differentiate functions of Grassmann 
variables. A partial derivative over 9 a is a linear operator 



d9 a Gn 



Qn 



which is defined by relations 



d 

and by the Leibniz's rule 
d 



d9 a 



i)t)a (0 b f(0))=5 ab f(9)-9 b — f(9). 



(7) 



It follows from this definition that partial derivatives 
must anticommute: 



d d 



d d 



39 a d9 b 89 b d9 a 



0. 



(8) 



Since a derivative g§-/(0) does not depend upon vari- 
able 9 a , it is sometimes convenient to think about dif- 
ferentiation as a linear operator mapping Q n into Q n -i 
(strictly speaking, after the differentiation one should re- 
name variables o +i, . . . , 9 n into 6 a , . . . , n -i)- Such an 
operator is called integration and is denoted as 

d9 a = —— : Q n — > Qn-i- 
We will also use a notation 

D9= I d9 n --- I d9 2 I d9 x . 



The order here is chosen such that J D9 6\ ■ ■ ■ 6 n = 1. It 
follows from Eq. [H] that 



so by applying the Leibniz's rule Eq. [7] one can use an- 
ticommuting version of integration by parts. 
Under a change of variables 



V a =J2 T *b0b, T G GL(n, C) 
6=1 



(9) 



operators of differentiation and integration change as fol- 
lows: 



2» 



" -in " 



6=1 



ba 89 b ' 



and 



Dr; = (dctT)" 1 D9. 



(10) 



(11) 



In the rest of this section we consider only even- 
dimensional Grassmann algebras Qm- Let us consider 
a quadratic form 

2n 

9 T M9 = M ab 9 a 9 b G Q 2n 

a, 6=1 

defined for any complex antisymmetric In x In matrix 
M and a bilinear form 

In 

9 T n = ^29 a n a G Q 4n - 

0=1 

Throughout the paper we shall extensively use these two 
formulas for Gaussian integrals: 



and 



D6»exp ( -9 1 M9 ) = i"Pf(M), (12) 



D0 exp ( rf 9 + -9 1 M9 ) = i n Pf (M) 



ex P ( ^2 r,TM 



(13) 
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In these formulas Pf(JV) is the Pfaffian of a complex an- 
tisymmetric matrix N defined as 'antisymmetrized prod- 
uct' A(Ni i2 N 3A ■ ■ ■ N 2n ,2n-i) that is 



Pi{N)= 2^ S sgn(a)iV CTl , CT2 ... 



Note that the integral in Eq. is an element of Q 2 



IV. GAUSSIAN STATES 

To any linear operator X G C2 n one can naturally as- 
sign a polynomial uj{X, 6) G Q 2n of 2n Grassmann vari- 
ables defined by 



uj(c p c q 



■ Of, I 



9 n 



1. 



This definition extends by linearity to an arbitrary X G 
C2n- We shall call o;(X, 6) a Grassmann representation 
of the operator X. It should be emphasized that uj is 
just an isomorphism of linear spaces. It has nothing to 
do with multiplication in the algebras C 2n and Q 2n - 

Consider as an example the projector d\a\ projecting 
onto a state "the first mode is empty" . In terms of gen- 
erators c one has 

aia\ = - (l+ icxc-i^ , 
Its Grassmann representation looks as 

co( ai al 8) = \{l + tdifc) = \ cxp (i0i0a). (14) 

Here and below all exponents are defined by their Taylor 
series. 

For any operators X, Y G C 2n one can compute a trace 
Tr(Xy) using the following simple formula 



Tr(XF) = (-2)" J DODp e° T,t u(X, 9)u(Y, p), 



(15) 



which can be verified by a direct inspection. 

Let V G G c be a canonical transformation implement- 
ing a rotation i? G SO(2n), see Eq. 0, and X G C 2n be 
an arbitrary operator. One can easily check that 



2)1 



;(VXV\8)=uj(X,i 1 ), Va = J2 R * 



fcCfc- 



(16) 



6=1 



Thus canonical transformations are equivalent to an or- 
thogonal change of a basis in the space of the Grassmann 
variables. We are now ready to define an important class 
of Gaussian states. 

Definition 1. A quantum state of n fermionic modes is 
Gaussian iff its density operator p G C 2n has a Gaussian 
Grassmann representation, that is 



■ exp 



M6 



for some 2n x 2n real antisymmetric matrix M. The 
matrix M is called a correlation matrix of p. 



Note that p itself generally can not be written as an 
exponent of an operator, since it might not have a full 
rank. By definition all Gaussian states are described by 
even polynomial in the Grassmann representation, see a 
remark [l4j . If p G C 2n is a Gaussian state, its correlation 
matrix can be found from 



-Tr (p[c a ,c b ]) = 



Tr (pic a c b ) for a ^ b, 
for a = b. 



All higher correlators are determined by the Wick for- 
mula, namely 

Tr(pFc ai c a2 ---c 2p )=Pi(M\ au ), (17) 



is 



with 1 < a\ < ■■■ < a 2p < 2n. Here M\ a 
the 2p x 2p submatrix of M with the indicated rows and 
columns. For example, 

Tr Goi 2 cic 2 c 3 c 4 ) = M 12 M U - M 13 M 24 + M U M 23 . 

It is known that any real 2n x 2n antisymmetric matrix 
M can be transformed by the adjoint action of SO(2n) 
into a block-diagonal form, with all blocks being 2x2 
antisymmetric matrices: 



Af = i?0 

3=1 



\ 3 
-A, 



R T , ReSO(2n). (18) 



The absolute values |Ai|, . . . , |A„| are referred to as 
Williamson eigenvalues of M. It follows from Eq. [S] that 
any Gaussian state p can be transformed by a canonical 
transformation into a product form 



1 - 

p(Ai, . . . , A n ) = — JJ(7 + i\jC 2j -i 



c 2j) 



(19) 



3=1 



Thus non-negativity of p is equivalent to having inequal- 
ities 



AjGt-1,1], 3 = 1,. 



,,n, 



or, in terms of operators, M T M < I. They are necessary 
and sufficient for a matrix M to be a correlation matrix 
of some Gaussian state. The state p is pure iff Xj = ±1, 
or M T M = I. Note that the density operator of the Fock 
vacuum is given by Eq. with all Xj = 1. 

Consider a projective measurement of an occupation 
number of some fermionic mode, say the first one, per- 
formed on a Gaussian state p with a correlation matrix 
M. Let us derive in few lines an explicit expression for 
the measurement outcomes probability distribution. The 
probability of the outcome "the first mode is empty" is 

Tr(paial) = (-2)" / L>^e 0TM w(aia|,6»V(p,/i), 



see Eq. |15j . It is convenient to introduce 2n x 2n matrix 
K with the only non-zero matrix elements K\ 2 = 1 and 
K 2 i = — 1, such that 

u{a 1 a\,e) = ^p(U T Ke\ 
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see Eq. |TI] . Using the rules Eq. |12I13| for taking Gaus- 
sian integrals, one gets: 

(-1)" 



Pf(M) / D0ei eT ^ K - M ~^ e 



(20) 
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= - Pf(M) Pf(i^ - M _1 ) 



If M is a singular matrix the last expression may be reg- 
ularized taking into account that Pf (L) 2 = det (L) for 
any antisymmetric matrix L. Thus we arrive to 



Tr(/9aia{; 



1 



det {MK - I), 



(21) 



where I is 2n x 2n unital matrix. A formula similar to 
Eq. |21j for an overlap between two fermionic Gaussian 
states has been known for many years, see 9]. In the rest 
of the paper we will write X(9) instead of u(X, 9) for the 
Grassmann representation of any operator X 6 C 2n . 



V. GAUSSIAN OPERATORS 

The main goal of this section is to prove that the FLO 
projectors map Gaussian states into Gaussian states. 
This fact has been already proved in the paper Q by 
demonstrating that each FLO projector is a limiting 
point for a converging sequence of generalized canonical 
transformations (implementing complex rotations from 
the group SO(2n, C)) and by using continuity arguments. 
We shall present an independent proof that relies on a re- 
markable symmetry of Gaussian states. 

Definition 2. 

(a) An operator X G C 2n with Ti(X) ^ is called Gaus- 
sian iff its Grassmann representation is 



X{6) = Cexp ( -9 T M9 



(22) 



for some complex number C and some 2n x 2n complex 
antisymmetric matrix M . We shall call M a correlation 
matrix of X . 

(b) An operator X G C 2n with Tr(X) = is called Gaus- 
sian iff 



X = lim X n 



(23) 



for some converging sequence of Gaussian operators 
X m G C 2n , Tr(X m ) ^ 0. 

In other words the set of Gaussian operators is a clo- 
sure of a manifold specified by Eq. [22]. For example, 
Gaussian states and FLO projectors (see Eq. \14\) are 
Gaussian operators. As we shall see, canonical transfor- 
mations are also Gaussian operators, but some of them 



have a zero trace and admit only a singular representa- 
tion as in Eq. EH (e.g. U = c x c 2 G G c , Tr (U) = 0). 
Roughly speaking, we have a regular representation for 
density operators and a singular one for operators de- 
scribing some transformations. (This is not a matter of 
concern, as these transformations also admit a regular 
representation in terms of Gaussian integrals, see Sec- 
tion EH) Note that all Gaussian operators are even. 

To treat all Gaussian operators on equal footing we 
shall give an equivalent definition which is less explicit 
but more operational. Consider two copies of the original 
system and an operator 



in 



(24) 



which belongs to the tensor product Ci n ® Ci n of two 
Clifford algebras 01 . One can easily verify that A is 
invariant under any canonical transformation V G G c 
applied simultaneously to both factors: 



[A, V ® V] = for any V G G c 



(25) 



Besides, a direct inspection shows that [A, p ® p] =0 for 
any 'product' Gaussian state as in Eq. |19|. Since an 
arbitrary Gaussian state can be converted to a 'product' 
one by a canonical transformation, it follows from Eq. |25j 
that [A, p ® p] — for any Gaussian state p. It turns out 
that the analogous commutation relation holds for all 
Gaussian operators and only for them. 

Theorem 1. 

(1) An operator X G Ci n is Gaussian iff X is even and 
satisfies 

[k,X®X] = 0. 

(2) Any Gaussian operator X has a form 



2k 



X{6) = C JJ /ia exp 



\a=l 



where p a = 53&=i T a },9b for some invertible complex ma- 
trix T . 

Let us first discuss several trivial consequences of the 
theorem. The part (1) implies that the set of Gaussian 
operators is closed under multiplication of operators: 

Corollary 1. IfX,Y G C 2n are Gaussian operators then 
XY G C 2n is also a Gaussian operator. 

Since the set of Gaussian operators is closed under tak- 
ing Hermitian conjugation, we infer that 

Corollary 2. If X is a Gaussian operator and p is a 
Gaussian state such that Xp ^ then a state 

p' = XpXy Tr(XpX^) 

is a Gaussian one. 




M ab p a p b (26) 
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This is quite important result since it tells that an 
arbitrary sequence of FLO operations can only produce 
a Gaussian state, provided that the initial state was a 
Gaussian one. To keep track of the state's evolution one 
just needs to compute a correlation matrix of p' given 
correlation matrices of p and X. We shall postpone a 
solution of this problem until Section IVT1 

In the rest of this section we prove Theorem ^ For 
convenience the proof is split into two lemmas. 

Lemma 1. For any Gaussian operator X £ Cm one has 



[A,X®X] = 0. 



(27) 



Proof: Applying the isomorphism of Section IIVI sepa- 
rately to both subsystems we get Cm ® C2n — Q2n ® Qim 
see [l3|. Our first goal is to describe the adjoint action 

A ad : 0->[A,0], 

which maps C2n <8 Cm into itself, in terms of Grassmann 
variables. Introduce a differential operator 

We claim that 

[c a ®c a ,Y ®Z}{9) = A a -Y(9)®Z(9) (28) 

for any operators Y, Z £ Cm having the same parity (i.e. 
for both Y and Z being odd or both being even). With- 
out loss of generality, both Y and Z are monomials in 
c's. In this case each of them either commutes or anti- 
commutes with c a . Consider two cases, (a) Both Y and 
Z contain c a , or both Y and Z do not contain c a . Then 
the commutator [c a ® c a , Y ® Z] is zero since both fac- 
tors yield the same sign (recall that Y and Z have the 
same parity) . The righthand side of Eq. [33] is also zero, 
since either 8 a or annihilates both Y and Z. (b) Y 
contains c while Z does not contain c a (or vice verse). 
In this case c a ® c a anticommutes with Y ® Z. Let us 
write Y = c a Y , where Y is a monomial which does not 
contain c a . We have: 

[ca ®c a ,Y ®Z] = 2{c a ® c a )(Y ® Z) = 2Y ® (c a Z). 

On the other hand, 

d 



d9 a 



Y®Z = 0, 



and 



— ®Q a -Y ®Z = Y ® 9 a Z. 

Substituting the last three formulas into Eq. [2H| we again 
get an equality. Summarizing, the adjoint action of A can 
be described by a differential operator 



2n 

A ad = 2^ (e a . 

a=l 



d d 



(29) 



Now we can easily conclude the proof. If A is a Gaus- 
sian operator with Tr(A) ^ 0, we can write 



X{6) = exp ( -9 T M6 



(the constant factor C is omitted) and thus 

/ 2n \ 

O 



w x = i[^M ab e b )x. 



\b=l 



Applying the differential representation of Eq. we get 

In 

A ad • X ® X = ( M ab + Mba) a ® b 

a,b=l 
■ X® X = 0. 

(30) 

Thus [A, A ® X] = 0. If A is a limiting point for a 
sequence of Gaussian operators {X m }, Tr(A m ) ^ 0, con- 
tinuity arguments show that [A, X ® X] =0 as well. □ 

In the next lemma we prove the reverse statement and 
as a byproduct derive an explicit expression Eq. j2Bl for 
traceless Gaussian operators. 

Lemma 2. Suppose X G Cm is an even operator such 
that [A, X ® X] — 0. Then X is a Gaussian operator. 

Proof: Let us first prove the lemma for the case Tr(A) ^ 
0. Making use of Eq. we can rewrite the equality 

[A, X ® X] = as 



o=l 



Denote C = 2~ n Tr (A) = A(0) and represent X{6) as 
X(9) = C ■ 1+ % — ^2 M ab a 6 b + higher order terms. 

a, 6=1 

Applying a differential operator 1 ® to both sides of 
Eq. pn] and making use of the Leibniz's rule Eq. [2] one 
gets: 



a— 1 x 

d 



(32) 



+ w b x ® x = °- 



Now let us put 9 = in the second factor (e.g. pick up a 
coefficient before 1). Taking into account the decompo- 
sition for X(9) given above one easily gets: 



_d_ 

Wb 



x = i^M ba e a x. 



(33) 
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This differential equation can be easily solved by making 
a canonical transformation which brings M to the block- 
diagonal form, see Eq. [18|. and using formulas Eq. |1U| . 
Taking into account "initial conditions" X(0) = C one 
gets 



X(0) = Cexp i^-9 T M9j, 

that is X is a Gaussian operator. 

Now let us prove the lemma for the general case. Con- 
sider a linear subspace jWi C Q 2n spanned by linear func- 
tions of Grassmann variables: 

2n 

Mi = {feg 2n ■■ f(e) = J2 a - e "}- 



Here the coefficients a a are complex numbers. Denote 
K, C A4i a subspace spanned by linear functions which 
annihilate X, i.e. 

JC = {fEM 1 : f(e)X(6) = 0}. 

(Since X is an even operator, any linear function / £ Mi 
commutes with X, so the left annihilation is equivalent 
to the right annihilation.) Let us perform a linear change 
of variables 



Ma — ^ T a b8b, 



(34) 



6=1 



with T being an invertible complex matrix chosen such 
that the first k variables fi span the subspace tC, i.e. 

K, = linear span \p\ , . . . , /Lt&] . 

From equalities [ijX = 0, j G [1, k] it follows that 

/ k \ 



(35) 



where X(fi) depends only upon Mfe+i! . . . , Hin- From 
Eq. we infer that 



In, q In q 

6=1 6=1 ' 



see Eq. [10]. It implies that the differential operator A a d 
is invariant under the change of variables: 



A 



ad = 2 21 Ma ® ^ h ^ <8> Ma 

Q=1 V "Ma "Ma 

Thus the function X(/i) satisfies an equation 



2d 



J2 i^^TT- + 7T- ®fJ-a)X®X = 0. (36) 
a=fc+i v 0fIa 0fIa 



Note that polynomials fi a X, a e [k+1, 2n] are linearly in- 
dependent, since otherwise we could extend the subspace 
JC. Then it follows from Eq. that the derivatives 
-^-X must be some linear combinations of polynomials 

fi a X, a £ [k + l,2n]. That is X(fi) obeys a differential 
equation similar to Eq. |33j and we conclude that 



2n 



X(/i)=Cexp|- ^ 

a,6=fc+l 



[l a M ah [l h 



for some complex (2n — k) x (2n — /c) matrix M and some 
complex number C. Since X was supposed to be an even 
operator we infer that k must be an even integer. Com- 
bining the formula above for X(^jl) and Eq. |35j we arrive 
to the desired representation Eq. for the operator X. 

If k = the lemma has been already proved. If k > 0, 
X is a limiting point for a sequence of regular Gaussian 
operators as in Eq. |221 since 

k 

nH a = lim exp (£/Lti/z 2 H htMfc-iMfc)- 

a=l 

We have proved that X is a Gaussian operator. □ 



VI. GAUSSIAN LINEAR MAPS 

By a linear map we shall mean a linear transformation 
sending operators to operators. A transformation of a 
quantum state under any physical operation (or under a 
conditioned physical operation such as a projector) can 
be described in terms of completely positive (CP) linear 
map. We have shown in the previous section (see Corol- 
lary 2 next to Theorem 1) that the adjoint action of any 
Gaussian operator is a CP map that has the following 
nice property: Gaussian operators are transformed into 
Gaussian operators. The main goal of this section is to 
introduce more general class of CP maps sharing this 
property and to describe explicitly their action on the 
level of correlation matrices. 

Definition 3. A linear map £ : Cm — > C-2n is Gaussian 
iff it admits an integral representation 

£(X)(6)=C [ exp[S{e,r ! )+ir 1 T ^]X(fi)DT 1 T)fi, (37) 



wht 



(38) 



for some complex In x In matrices A, B, D, and some 
complex number C . 

The function S(9,rj) will be called an action for the 
map £ . Since the matrices A and D enter into Eq. |38j 
as 9 T AO and rj T Drj, we can assume them to be antisym- 
metric. Let us emphasize that no restrictions are put 
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here on the input operator X\ it may be non-Gaussian 
and need not to have a definite parity. The double inte- 
gration in Eq. |37| may be thought of as a trace over one 
copy of the system, see Eq. |T5] . 

As a simple example let us consider the identical map 
£i(X) = X. One can easily verify that its integral rep- 
resentation is as follows: 



p 2 71 

X{6) = j n(0a-A.)*(M)D/i. 



Taking into account that 



2n „ 

l[(9 a - ,i a ) = J exp [i(6 - nfrj\ D v , (39) 



we finally get: 

X(9) = J exp [i6 T i-i + iri T v\X{p)T)riT>ii. 

Thus Zi is a Gaussian map with an action Si(0,r)) = 
i8 T r], i.e. A = D = and B = I. 

A canonical transformation implementing a rotation 
R G SO(2n) of the generators c, see Eq. [5], is described 
by an action Sn(9,rj) — i8 T R T r], that is A = D = 
and B — R T . Indeed, according to the previous para- 
graph, it specifies a linear map £r such that £r(X)(6) — 
£ I (X)(R9) = X(R9). This orthogonal change of vari- 
ables is exactly the effect of canonical transformations in 
the Grassmann representation, see Eq. 

Now we shall establish several useful properties of 
Gaussian maps. Recall that an operator X e C 2 n is called 
even (odd) if it is a linear combination of only even (odd) 
monomials in the generators c. 

Definition 4. A linear map £ : Ci n — > Ci n is called 
parity preserving iff it maps even (odd) operators into 
even (odd) operators. 

Lemma 3. Any Gaussian map is parity preserving. 

Proof: For any operator X £ C2« define an operator X £ 
C-in such that X(6) = X(—9). Obviously, X is an even 
(odd) operator iff X = X (X = —X). Let us make a 
change of variables /i — » — /i and rj — > —77 in the integral 
Eq. |2Z1- Taking into account that S(0,-rj) = S(-9,rj) 
and Eq. [TTJ saying that D( -fi) = D(/i), D( -rj) = T>(rj) 
we arrive to 



£(X)=£(X). 



(40) 



Thus if X = ±X then £(X) = ±£(X). The lemma is 
proved. □ 

Lemma 4. Gaussian maps transform Gaussian opera- 
tors into Gaussian operators. 

Proof: Instead of proving this statement by direct appli- 
cation of Eqs. |12I13| for Gaussian integrals (which is not 



so simple due to possible singularities), we will make use 
of Theorem^ (part 1) and the differential representation 
Eq. f° r the adjoint action of the operator A. Since 
Gaussian maps are parity preserving and Gaussian oper- 
ators are even, one suffices to prove that 



A ad • £{X) £(X) = 



(41) 



for any Gaussian operator X and any Gaussian map £ . 
The proof involves three main ingredients: (i) an identity 

E£i ® A + A ® *«) ' exp (ie T B V f 2 



Ea=l (Va 



rj a ] • exp [i9 T Br]) 



which can be easily verified; (ii) the Leibniz's rule 



wj fg) = (§9z) 9+f (w a 



(42) 



(43) 



which is valid for even polynomial / and arbitrary g; (iii) 
the integration by parts formula, i.e. 

This formula is valid if one of the functions / and g is 
even. Here the sign '+' stands for even g and the sign 
'-' stands for even /. (If both functions are even then 
the integrals on the right and on the left are both equal 
zero, so one can choose an arbitrary sign.) Introduce an 
auxiliary polynomial 

f(r]) = J exp [iif^i] X{ji) D/i . 

By the same arguments as in the proof of Lemma [21 one 
can show that / is even. Let us apply A a( j to a polynomial 

P x = (I exp [S(9, 77)1/(77) Dryj e Q 2n ® Qm- 

We can move A a d from the left to the right applying sub- 
sequently (a) The Leibnitz's rule Eq. 021; (b) Lemma ^ 
(c) The identity Eq. 02; (d) The integration by parts 
formula Eq. 031; (f) The Leibnitz's rule; (g) LemmaC] 
After these steps one gets 

A, d Px = -[ exp[S(9,f 1 )f 2 A ad f(ij)'» 2 Di 1 ®T)f ] . 

Let us keep moving A a( j to the right by applying subse- 
quently the identity Eq. |42j and the integration by parts 
formula Eq. 031- After these moves one gets: 

A ad /(r^ 2 = - / exp [777 T '/jf 2 A ad A( M )^ 2 D/i ® D/a . 

We arrive to 

A ad ■ £{X) «) £{X) =£®£ (A ad ■ X ® X) = 0. 
The lemma is proved. □ 
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Consider a Gaussian operator X which can be de- 
scribed by a correlation matrix M as in Eq. |22j and 
a Gaussian map £ as in Eqs. |37I38| . Lemma 0] implies 
that £{X) is a Gaussian operator. Applying the Gaus- 
sian integration rule Eq. ^3] one can show that £ (X) has 
a correlation matrix 

£{M) = B (M _1 + D) 1 B T + A (45) 
= B (I + MD)' 1 MB T + A, 

while a pre-exponential factor of the operator £{X) can 
be found from an identity 

Tr (£(X)) = C(-l) n Pf(M) Pf(M" 1 + D) Tr {X). (46) 

The value of Tr (£(X)) can be found up to a factor ±1 
using a regularized version of Eq. : 

Tr {£ {X)f = C 2 det {I + MD) Tr (A) 2 . (47) 

The expression Eq. 03] does not make sense if (I + MD) 
is a singular operator, see a comment [l^. If this is the 
case, the Gaussian operator £(X) admits only a singular 
representation, as in Eq. [25] . 

Suppose that Tr (£(X)) = Tr (X) for any operator X, 
that is the map £ is trace preserving (TP). Then it fol- 
lows from Eq. g?| that a function f(M) = det (I + MD) 
does not depend on M, as long as M is an antisymmet- 
ric matrix. From the equality f{D^) = f(-D') we get 
det(/ + D t £)) = dct(/-D t D) which is possible only 
for D = 0. Then it follows from Eq. gH] that £ is TP 
iff D = and C = 1. A transformation of correlation 
matrices for TP maps looks particularly simple: 

£{M) = BMB T + A. 

A TP map £ preserving an identity operator, £{I) = /, 
is called bistochastic. Since the identity operator is a 
Gaussian one with a zero correlation matrix, we conclude 
that a Gaussian map £ is bistochastic iff A = D = 0, 
C = 1. A transformation of correlation matrices for bis- 
tochastic maps is just 

£{M) = BMB T . 

VII. COMPLETELY POSITIVE GAUSSIAN 
MAPS 

Our next goal is to work out conditions under which 
a Gaussian map is completely positive and thus may de- 
scribe a physical transformation of states (may be a con- 
ditional transformation such as a measurement with a 
postselection). By definition, a map £ : C 2n — > C 2n is 
completely positive (CP) iff a map 

£®Id : C 2n ® C 2n -> C 2n ® C 2 „ 

is positive. Here Id stands for the identity map. 



Theorem 2. A Gaussian map £ specified by Eq. f5 l 7l5 l #j ? 

is completely positive iff (i) C > and (ii) a matrix 

is real and satisfies M T M < I . 

Before proving the theorem we shall introduce a 
fermionic tensor product of linear maps <g>f and a no- 
tion of a dual state. Let Z £ C^ n be a monomial in the 
generators c. One can always represent it as Z — Z' ' Z" ', 
where Z' is a monomial in c\ , . . . , c 2n only and Z" is a 
monomial in c 2n +i, . . . , C4„ only. 

Definition 5. Let £\,£ 2 : C 2n C 2n be linear maps. 
Define a map £\®f£ 2 '■ C^ n — > C^ n on monomials Z 
'/"/" as 

(£i ®f£ 2 ){Z) = £ 1 {Z')£ 2 {Z"), 
and extend it using linearity to the whole algebra Ci n . 

In other words, the map £\ ®f£ 2 acts on the first half 
of the modes as £ i and on the second half of the modes 
as £ 2 . 

Lemma 5. Let £ be a parity preserving linear map. 
Then £ ® Id is positive iff £ ®fld is positive. 

Proof: Define a Jordan- Wigner transformation 
J : Ci n — > C 2n C 2n 

according to 

J(c fc ) = c k ®P for 1 < k < 2n, (48) 
J(cfc) = / <8) Cfc for 2n < k < An. 

Here P is the parity operator: 

P = i n c 1 c 2 ---c 2n , P 2 = I, P^ = P. 

The definition Eq. 0Sj extends to arbitrary operators 
by linearity and multiplicativity, J(XY) = J(X)J(Y), 
since J preserves commutation relations between the 
generators. Thus J is an isomorphism of algebras C^ n 
and C 2n (£> C 2n . Since J preserves Hermitian conjuga- 
tion rules, J(X^) = J(Xy, it can be represented as 
J(X) = WXW', where W is a unitary operator on the 
space of states of 2n fermionic modes (which is isomor- 
phic to the space of states of two copies of n fermions). 
For any monomial Z = Z' Z" £ C^ n as above one has 

J(Z'Z") = Z' ® P e Z", 

where e = for even Z' and e = 1 for odd Z'. Thus 

((£ ® Id) o J) (Z) = £{Z') ® P e Z". 

Besides, according to Definition|S]we have (£®fId){Z) = 
£(Z')Z" . Since £ is parity preserving, £(Z') is a linear 
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combination of monomials having parity e. Using linear- 
ity of J we get 

(J o (£® f Id)) (Z) = £{Z') ® P'Z". 

We conclude that 

(£®Id)oJ = J o(£® f Id). 

Since J is a unitary isomorphism of algebras, it follows 
that positivity of £ ® Id is equivalent to positivity of 
£® f Id. □ 

To check positivity of the map £ <£> / Id we shall make 
use of Jamiolkowski duality between maps and states, 
see ^JEI; slightly adapted to fermions. 

Definition 6. Let £ : Cm — > C 2n be a parity preserving 
linear map. An operator pg G C^ n dual to the map £ is 
defined " 



as 



P£ = (s® f id)( PI ), 



whe 



Pi 



^ 2n 

2^ n ( 



I + ic a C2n+a ) G C 



a=l 



Note that pi is a pure Gaussian state. It can be 
regarded as a 'maximally entangled' state between the 
modes t\, . . . , c 2n and the modes c 2n +i, . . . , c\ n . A dual 
operator completely specifies a linear map. To check that 
let us use the isomorphism C± n = Gin, see Section llVI 
Assign Grassmann variables 6\ , . . . , 9 2n to the genera- 
tors Ci , . . . , c 2n and variables 771 , . . . , rj 2n to the generators 
C2n+i, • • ■ , C4n- Then one has 



Pl( 



It follows that ps{6, rf) is just a generating function for 
the map £ , so by taking derivatives over r\ one can find 
an image under £ of any monomial in c generators. 

A Gaussian map defined by Eqs. |37I38| has a dual 
operator 



- J exp [S(0, C) + *C T M + *p T v] D£i DC . 



The integral over p gives the integral kernel for the iden- 
tity map, see Eq. The integration over £ yields 



PS(^V) = §^^P[S(0,V)]- 



(49) 



Thus Gaussian maps have Gaussian dual operators. Now 
we are ready to prove the main theorem. 

Proof of Theorem\^ (a) Necessity: If £ is a CP map, 
Lemma |5] implies that ps must be a non-negative op- 
erator. Since Trp,r = C, we conclude that C > 0. Since 
ps is a self- adjoint operator, all matrix elements of M 



must be real. Thus pg is proportional to a density op- 
erator of some Gaussian state, see Definition ^ that is 
M coincides with a correlation matrix of some Gaussian 
state. We have already shown in Section llVl that its non- 
negativity is equivalent to a constraint M T M < I on a 
correlation matrix. 

(b) Sufficiency: The conditions of the theorem are equiv- 
alent to non-negativity of pe- To prove complete positiv- 
ity of £ let us show that it admits a Kraus representation, 
that is £{X) = J2 a AaXA^ for some operators A a G C 2n - 
The proof is a trivial adaptation of standard "bosonic" 
arguments (see |ll| and references therein) to the case of 
fermions. 

A pure state |i>) such that pi = | < I ) )( < I ) | (unique up to 
a phase) satisfies 



C 2 «+a|$) = -ic a |$), 



1, 



,2n. 



It means that an arbitrary pure state of 2n fermionic 
modes can be written as |\&) = for some operator 

A G which involves only the modes c\ , . . . , c 2n from 
the first half. Consider an arbitrary pure state decompo- 
sition of ps ■ 



Ps 



E i*«x*«i- 



Representing |\& a ) = A a \&) we get 



ps 



\A aPl Al 



Here all A a include only the modes ci, . . . , c 2n . Since a 
dual state completely specifies a map, we conclude that 
£{X) — J2 a AaXA^ for any operator X G C 2n - Com- 
plete positivity of £ follows from existence of a Kraus 
representation, see ^2|- ^ 

As a simple application of the theorem, let us char- 
acterize bistochastic CP maps. As we already know, 
A = D = and C = 1 for any bistochastic map, see 
the last paragraph of Section I VII Then the complete 
positivity is equivalent to the matrix B being real and 
an inequality B T B < I. Thus the matrix B can be rep- 
resented as B = RiBR r , where Ri,R r G SO(2n) are 
rotations and B is a diagonal matrix, 

B = diag[Bi,B 2 , • • - ,-B 2n ], 

such that \B a \ < 1 for all a. The rotations Ri and R r 
can be undone by composing the map £ with appropriate 
unitary canonical transformations on the left and on the 
right, see Eq. [5]. We conclude that 

£ = £ 1 o £ o £ r , 

where £\ , £ r are unitary maps and £ has an integral rep- 
resentation 



£(X)(9) 



exp 



: - E B a® alia + VaPa 



X(p) DjjDp. 
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Evaluation of this integral yields the following action of 
£ on the monomials in c generators: 

£(Cai ' ' ' ^afc ) — B ai • • • B ak C ai • ' ' C ak j 

where k = 1, . . . , 2n. In other words, £ is a composition 
of elementary maps 



\BJ < 1, 



1. 



,2n. 



which all pairwisc commute. It is a direct analogue of a 
product map for qubits or bosonic modes. 

One can ask whether the set of Gaussian maps is closed 
under a composition of linear maps. As we shall see now, 
in many important cases the answer is affirmative. 

Lemma 6. Suppose £± and £ are trace preserving Gaus- 
sian maps. Then £ 2 °£i is a trace preserving Gaussian 
map. 

Proof: Denote £ = £% ° £% • It follows from Lemma 
that £ is a parity preserving map. Thus one can define 
a state ps dual to the map £ . Lemma 0] implies that 
ps is a Gaussian operator. Since Tr ps = 1 we conclude 
that ps admits a regular exponential representation as in 
Eq. and thus a map £ admits an integral represen- 
tation Eq. |37j . that is £ is a Gaussian map. □ 

Lemma 7. Suppose £\ and £2 are completely positive 
Gaussian maps. Then £2 ° £1 is a completely positive 
Gaussian map. 

Proof: The first part of the proof copies the one of the 
previous lemma. The dual state ps is a non-negative 
operator and thus Tr ps > unless ps = 0. In the latter 
case £ = 0, since a dual state completely specifies a linear 
map. The map sending all states to zero is obviously a 
Gaussian one. In the former case one can apply the same 
proof as above. □ 



VIII. FLO OPERATIONS AS GAUSSIAN MAPS 

In this section we put forward an integral representa- 
tion for a single-mode occupation number measurement. 
It will allow to find explicitly how do the FLO measure- 
ments change a correlation matrix of the measured state 
(for each of the two outcomes) . 

The measurement of the j-th mode is described by 
orthogonal projectors ajcij and djOj which correspond 
to the outcomes "the mode is occupied" and "the mode 
is empty" respectively. Let us introduce linear maps 
£j t o,£j t i : C2n - > C2n describing the adjoint action of 
the two projectors: 



(50) 



The index j runs from 1 to n, where n is the total number 
of fermionic modes. Our first goal is to find operators 



dual to these maps, see Definitional 

Pj.o = {£j,o ®fld)(pi) = (oj-ot) pi (ajot), 

and 

Pj.i = (£j,i ®fld)(pi) = (otoj-) pi (a]a 3 ). 

Here the operators a], j — 1, . . . ,n are considered as 
elements of the Clifford algebra C4„, since pi £ Ci n . After 
a simple algebra one gets: 

1 

Pjfi 



2 2n+T \ T + i£ 2j-lC2jJ [I - iC2n+2j-lC2n+2j 
■ Y[ (j + ic a C 2n +a 



Here the product \\ is taken over all a in the range 
1, . . . , 2n excluding a = 2j — 1 and a — 2j. Analogously 
one gets 

1 



2'2n+l 



(l - iC2j-lC 2 j) (j + iC2n+2j-lC2n+2j) 



n'(' 



+ ic a C2n+a 



Let us assign Grassmann variable 81, ... , 82n to the gen- 
erators ci , . . . , C2n and variables rji , ... , r\2 n to the gener- 
ators C2 n +i, . . . , Ci n - Then the states pj >€ , e £ {0, 1} have 
the following Grassmann representation: 

1 



where 



+ i^ OaVa- 



Here the sum ^ is taken over all a in the range 1, . . . , 2n 
excluding a = 2j — 1 and a = 2j. Since a dual state com- 
pletely specifies a linear map, it follows from Eq. 05] that 
£j^ is a Gaussian map with an integral representation 

£j, e (X)(e) = i y exp [S j>e (0,v) +iv T n]X(n)T>r 1 T) t i. 

If A is a Gaussian state with a given correlation matrix 
then £j >€ (X) is a Gaussian state whose correlation matrix 
is given by Eq. |45|. 

Consider as an example a measurement of the first 
mode and the outcome "the mode is empty" . Let p £ C211 
be the input Gaussian state with a correlation matrix M, 

p(0) = ±-ex V (^8 T Me 

Computing Gaussian integrals in Eq. |51j one gets: 



p( ai a\) (8) = p-exp 



2n 



18182 + — ^ L a b9 a t 



a, b— 3 



(I - MK)~ X M. 
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Here K is an antisymmetric 2n x In matrix with the only 
non-zero matrix elements if 12 = —K21 = 1, that is 

K pq = S p lSq2 — 8p26 q l. 

The coefficient p turns out to be 

Po = ipf(Af)Pf(X-Af- 1 ), 

which agrees with the earlier derived formula Eq. |21j for 
the probability to observe the mode in the empty state. 

In general, if a state p to be measured is Gaussian with 
a given correlation matrix M, one can use the integral 
representation Eq. [SI] and the equation Eq. [15] to find 
a correlation matrix £j^ n (M) of the post-measurement 
state provided that the outcome was n. The probabilities 
of the two outcomes are given by Eqs. |46l4Tj . 

Evolution of a correlation matrix under unitary ele- 
ments of FLO (canonical transformations) is given by 
Eq. |16j . Thus the technique developed in the paper al- 
lows to simulate FLO on a classical probabilistic com- 
puter and both aspects of the simulation mentioned in 
the introduction can be easily addressed. 



IX. CONCLUSION 

Notions of a Gaussian state and a Gaussian linear map 
are generalized to the case of anticommuting (Grass- 
mann) variables. Conditions under which a Gaussian 
map is trace preserving and (or) completely positive are 
explicitly described. This formalism allows to develop the 
Lagrangian representation for FLO, which covers both 
unitary operations and projectors describing single-mode 
measurements. Using the Lagrangian representation we 
have reduced classical simulation of FLO to computation 
of Gaussian integrals over Grassmann variables. Explicit 
formulas describing evolution of a quantum state under 
FLO operations have been put forward. 
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